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Nonlinear interaction of a low density electron beam with a uniform plasma is studied using 
two-dimensional particle-in-cell (PIC) simulations. We focus on formation of coherent phase space 
structures in the case, when a wide two-dimensional wave spectrum is driven unstable, and we also 
study how nonlinear evolution of these structures is affected by the external magnetic field. In the 
case of isotropic plasma, nonlinear buildup of filamentation modes due to the combined effects of two- 
stream and oblique instabilities is found to exist and growth mechanisms of secondary instabilities 
destroying the BGK-type nonlinear wave are identified. In the weak magnetic field, the energy of 
beam-excited plasma waves at the nonlinear stage of beam-plasma interaction goes predominantly 
to the short- wavelength upper- hybrid waves propagating parallel to the magnetic field, whereas in 
the strong magnetic field the spectral energy is transferred to the electrostatic whistlers with oblique 
propagation. 
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INTRODUCTION 



Investigations of various aspects of the beam-plasma 
instability are important for our understanding of physi- 
cal processes occurring in space and laboratory plasmas. 
In recent years, regimes of collective beam-plasma inter- 
action typical to the fast ignition scheme or space and 
ionospheric phenomena have attracted particular inter- 
est. That is why most of recent theoretical and nu- 
merical studies deals with strong instabilities excited by 
counterstreaming electron beams of comparable densi- 
ties. Indeed, in the case of isotropic plasma, results of 
linear theory relevant for such regimes have been re- 
visited and dominance of different unstable modes in 
parameter space have been determined J] H]. Non- 
linear effects such as beam trapping and coupling of 
the Weibel-filamentation instability with the two-stream 
instability have been studied in the framework of the 
initial-condition problem using both two- [SHS] and three- 
dimensional PIC simulations. As to the case of mag- 
netized plasma, the particular attention has been focused 
on formation and stability of nonlinear phase space struc- 
tures (electron holes and tubes) arising during saturation 
of the two-stream instability j9HT4]. 

The goal of this paper is to investigate in details the 
nonlinear evolution scenario of a beam-plasma system in 
the case of low density beams <C Up and to study how 
such a scenario is modified by the external magnetic field. 
Qualitatively, the overall picture of weak beam-plasma 
interaction resembles that for a strong beam, but, in our 
opinion, more attention should be given to the study of 
consecutive energy transfer from one plasma mode to an- 
other. Our interest to this problem is motivated by the 
fact that wave-wave interaction due to the beam non- 
linearity that is responsible for the energy exchange be- 
tween plasma modes in a uniform beam-plasma system 
does also play an important role in the realistic problem 



of beam injection through a plasma boundary. In one- 
dimensional simulations [15 we have found the regime 
of beam relaxation, which is characterized by the reg- 
ular nonlinear beam dynamics in electric fields of res- 
onant waves even in a quasistationary turbulent state. 
The beam nonlinearity in this regime results in growth 
of coherent wave packets instead of waves with chaotic 
phases and leads to saturation of the pumping power. 
The key argument against feasibility of such a regime is 
a one-dimensional character of our simulations. Thus, it 
is interesting to study whether nonlinear beam dynamics 
plays an essential role in multidimensional problems and 
whether multimode character of the beam-plasma insta- 
bility is able to destroy correlation effects. For this pur- 
pose, we simulate two-dimensional beam-plasma inter- 
action in the framework of the periodic boundary prob- 
lem, in which different resonant modes interact with each 
other for a long time due to the beam nonlinearity. 

In Sec. |ll|we choose physical parameters of the beam- 
plasma system and numerical parameters of our PIC 
model. In Sec. |III| we study the linear stage of the beam- 
plasma instability and compare simulation results with 



theoretical predictions. Sec. |IV|presents simulation stud- 
ies for the nonlinear stage of beam-plasma interaction. 
Here, we identify the main nonlinear processes respon- 
sible for the energy exchange between different plasma 
modes and study the magnetic field effects. In final Sec. 
Mwe discuss our main results. 



II. SIMULATION PARAMETERS 

Let us consider relaxation of an electron beam with 
the density ni^/np = 0.002 and the temperature T5 = 10 
eV propagating with the velocity v^/c = 0.382 (c is the 
speed of light) in a uniform plasma with the initial elec- 
tron temperature Tg = 60 eV. This set of parameters is of 
a great importance for laboratory experiments on turbu- 
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Figure 1: (a)-(c) Theoretical growth-rates of unstable modes in the beam-plasma system for different values of external magnetic 
field, (d)-(f) Time-evolution of Ex and Bz Fourier harmonics at the linear stage of beam-plasma interaction in PIC simulations. 



lent plasma heating in mirror traps [TG'.TT^. Since plasma 
in such traps is confined by the external magnetic field 
B that can be characterized by the ratio of the cyclotron 
ujc = eB / {rriec) to the plasma ujp = {Ane'^np/mey^'^ fre- 
quency of plasma electrons Q = ujc/^p^ it is interesting 
to consider different cases Q = 0, Q = 0.5 and Q = 2 {e 
and rrie are the charge and the mass of electrons). 

To simulate nonlinear evolution of the beam-plasma 
system, we use the standard two-dimensional electromag- 
netic PIC code with periodic boundary conditions for 
particles and fields. Initially, the charge and current of 
the electron beam propagating along x are compensated 
by charges and currents of plasma electrons. We choose 
the simulation box x Ly = 360Ax x SSAAy and grid 
sizes Ax = Ay = 0.02c/ ujp in such a way as to fulfill 
the Cherenkov condition for the modes with wavenum- 
bers kx = Qn/Lx = oOp/vi) and resolve the Debye length. 
For the time step we get r = O.Olcj"^ and for the rea- 
sonable level of noise we use 256 computation particles 
in each cell. In addition to beam and plasma electrons, 
we take into account the dynamics of ions with the mass 
rrij = 1836mp. 



III. LINEAR ANALYSIS 

Let us find out what plasma modes should be driven 
unstable in the beam-plasma system according to the lin- 
ear theory and how their growth-rates depend on the 
external magnetic field . Due to low beam and plasma 



temperatures the linear analysis can be restricted by the 
fluid approximation. In this limit, the eigenfrequencies 
of unstable plasma oscillations can be found from the 
equation 
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where eai3 is the dielectric tensor with the components: 
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Here, 7 is the relativistic factor of the beam. In our axis 
the beam velocity and the magnetic field are directed 
along X and the wavevector has the components k = 

In the spatially periodic system x Ly possible 
wave numbers are restricted by the following discrete sets: 
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hence, the position of each plasma mode in the plane 
{kx, ky) can be marked by two integers n-m. In our case, 
the unstable spectrum contains only modes with numbers 
n = 0, 1, 2, 3. Numerical solutions of the dispersion rela- 
tion for the growth-rates of these modes as functions of 
the transverse wavenumber ky are shown in Fig. [TJa)-(c). 
It is seen that the transition from the case of isotropic 
plasma with dominant oblique instabilities to the case 
of strong magnetic field is accompanied by the essential 
reduction in the growth-rates of oblique modes and com- 
plete stabilization of the filamentation instability. 

Let us see how accurately theoretical predictions are 
reproduced in PIC simulations. Fig. [l](d)-(f) show that 
Fourier harmonics of Ex and Bz fields does really demon- 
strate exponential buildup, but not all of them grow with 
theoretical growth-rates. In numerical simulations, dom- 
inant modes appear to have an essential impact on slower 
instabilities with the same transverse wavenumber. In- 
deed, as one can see from Fig. [l] (d), in the case of 
isotropic plasma, growth-rates of dominant Langmuir 3- 
m modes agree well with theoretical predictions up to 
large m thus confirming fiow pattern of the observed in- 
stability, whereas instabilities of beam modes 1-m and 
2-m are found to be suppressed. As to the pure trans- 
verse filamentary perturbations, it is shown, using the 
0-3 mode as an example, that barely visible exponential 
growth of these modes with linear growth-rates switches 
rapidly to the nonlinear regime with effective rates ex- 
ceeding the maximum growth-rate of the linear theory. 
Nonlinear nature of this regime is confirmed by the fact 
that by the moment Upt = 75 fiuctuations of beam den- 
sity become comparable in magnitude with the unper- 
turbed value. 

The weak magnetic field Q = 0.5 reduces growth-rates 
of oblique upper-hybrid modes 3-m in a good agreement 
with the linear theory and keeps beam modes 1-m and 
2-m stable. In the strong magnetic field the spectrum of 
unstable plasma waves narrows significantly. It results in 
the situation when, in addition to the 3-0 and 3-1 modes, 
we observe the instability of the 2-2 mode. In contrast 
to previous cases, this mode now falls in the ky range, 
inside which it becomes dominant. 



IV. NONLINEAR STAGE 

Before proceeding to the detailed study of nonlinear 
evolution scenarios of the beam-plasma system in three 
different cases Q = 0^ Q = 0.5 and ^7 = 2, let us briefiy 



describe the main stages of beam-plasma interaction by 
considering time evolution of dominant Fourier harmon- 
ics of the electric field Ex presented in Fig. |2j In the 
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Figure 2: (Color online) Time evolution of Ex Fourier har- 
monics at the nonlinear stage of beam plasma interaction in 
different magnetic fields 11 = (a), Vi — 0.5 (b) and Vi — 2 
(c). 

case 17 = 0, at the stage of beam trapping, the highest 
level of saturation is achieved by modes 3-0 and 3-1. The 
resulting nonlinear wave, however, turns out to be unsta- 
ble against perturbations with shorter longitudinal wave- 
lengths (n = 4 and n = 5). Among these primarily stable 
plasma oscillations, the main role is played by the Lang- 
muir waves with small propagation angles (m = 0,1). 
Fig. [2] (a) shows that the time-averaged growth of these 
modes is found to be exponential and saturation levels 
of secondary instabilities appear to be so high that from 
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Figure 3: (Color online) Evolution of the beam density map (top row) and beam phase spaces (x,pcc) (middle row) and {x,py) 
(bottom row) in the case of isotropic plasma. Both phase spaces are obtained by averaging on y inside the region restricted by 
horizontal lines in the beam density map. Lines in the phase space {x,px) show px corresponding to phase velocities of different 
modes: 3-0 (solid), 4-0 (dashed) and 5-0 (dash-dotted). 



some time the spectrum of field is dominated by 4- 
and 4-1 modes. Further damping of these modes is 
accompanied by the modulation instability resulting in 
absorption of their wave energy by plasma electrons. 

In the weak magnetic field Q = 0.5, the picture of 
nonlinear evolution does not change drastically. Upon 
saturation of primary unstable upper-hybrid modes 3-0 
and 3-1, we observe the buildup of secondary instabili- 
ties with almost the same averaged growth-rates as in the 
previous case (Fig. [2](b)). At the saturation state, how- 
ever, secondary upper-hybrid modes with n = 4 do not 
play the dominant role. Significant changes in the picture 
of spectral energy transfer are found to occur only in the 
strong magnetic field Q = 2. In this case, the buildup of 
initially stable oblique modes 3-m with m > 1 becomes 
more efficient than the growth of longitudinal oscillations 
with n = 4 and n = 5 (Fig. [5] (c)). 



A. Isotropic plasma = 

To study nonlinear beam-plasma interaction in more 
details, we present time evolution of the beam density 
map ni){x^y) as well as snapshots of beam phase spaces 
(x,p^) and {x,py) (Fig. |3|. 

Beam trapping in oblique directions is the first non- 



linear consequence of the prevailing growth of almost 
electrostatic oblique modes 3-m. Indeed, in the time 
interval uOpt = 75.1 ^ 100.5, Fig. |3] demonstrates the 
process of beam density modulation under the most un- 
stable oblique modes, which is followed by the formation 
of characteristic vortex structures in the phase space of 
beam electrons. As we have mentioned, at the beginning 
of beam modulation we observe that the filamentation 
instability switches to the nonlinear regime. Let us dis- 
cuss the possible mechanism of nonlinear excitation of 
filamentary perturbations. 

Modulation of the beam density under unstable modes 
3-m and 3-(-m) drives density perturbations with the fol- 
lowing structure: 

Sub oc e^^^^(^/^^ - ^) + ^s-mt cos{27rym/Ly). 

Here, we neglect the difference between the eigenfre- 
quency of the mode 3-m and the plasma frequency ujp. 
Since at the same time the mode 3-0 perturbs the parallel 
velocity of the beam 

the parallel beam current contains the nonlinear term 
djx (X 6nbSvb (X e*^^3-m + ^3-o)t cos{27rym / Ly) , 
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which generates the magnetic field Bz with the same spa- 
tial structure. Thus, at the early stage, when density and 
velocity fluctuations are small enough, fllamentary per- 
turbations 0-m due to the combined effects of two-stream 
and oblique instabilities should grow with the approxi- 
mate rate 2r3_o. In our simulations the corresponding 
growth-rate is slightly lower. Apparently, it is explained 
by the fact that growing beam perturbations at some 
time cease to be weakly nonlinear. 

As it is seen from Fig. |3) after the trapping stage the 
beam breaks up into bunches localized in space in both 
parallel and perpendicular directions. One can also see 
that these bunches are arranged one behind the other 
thus forming current layers along the magnetic fleld. 
Since the beam acquires the large spread in transverse 
momentum under oblique modes, different current lay- 
ers exchange beam particles and merge into one by the 
moment ujpt = 252.5. At the same time, in phase spaces 
{x^Px) and (x^py) we observe the formation of nonlinear 
BGK-wave, the amplitude of which, in contrast to the 
one-dimensional problem, is localized in the transverse 
direction. 

Let us investigate stability of the resulting nonlinear 
wave. As we have mentioned, at this stage, the beam- 
plasma system is unstable against short-wavelength al- 
most longitudinal Langmuir oscillations with n = 4 and 
n = 5. The exponential growth of secondary instabilities 
indicates that the initial stage of these instabilities can 
be described by the theory linearized in the amplitudes 
of unstable perturbations. In this case, however, we can- 
not consider a uniform beam-plasma system even with 
the increased momentum spread of the beam as an un- 
perturbed state. Excitation of secondary plasma modes 
differs markedly from Cherenkov mechanism. Fig. [2] and 
[3] show that the instability of the mode 5-0, for exam- 
ple, starts to grow when there are no any particles in 
resonance with this mode (the momentum correspond- 
ing to the phase velocity of this mode is shown in the 
phase space {x^Px) as the dash-dotted line). The stabil- 
ity analysis of a one-dimensional (uniform along y) non- 
linear BGK-wave as an exact stationary solution of the 
Vlasov-Poisson system seems to be a more adequate theo- 
retical approach in our case. According to the theory pro- 
posed in Ref. [18, W, efficient energy exchange between 
a plasma mode with the frequency uj and the wavenum- 
ber k and a nonlinear BGK-wave with the phase velocity 
Vq and the bounce-frequency occurs if the following 
resonance condition is fulfllled: u = kvo + Nuij. Let 
us compute growth-rates of sideband modes 4-0 and 5-0 
driven unstable due to resonances N = —1 and A/" = — 3, 
respectively. 

For this purpose, we simplify the general theory of Ref. 
[19] using the following assumptions. Most of beam parti- 
cles is concentrated near the bottom of the potential well 
with the parabolic proflle V{x) = —e(po{l — k^x'^ /2) and 
their distribution function f(W) is constant inside the en- 
ergy range — e(/9o <W < —{).^eipQ = Wq ((/?o denotes the 
amplitude of electrostatic potential in the BGK-wave). 



In this case, all beam particles participate in bounce os- 
cillations with the same frequency = ko{e(po/m)^^'^ 
and their maximum deviation length from the center of 
the potential well is determined by 

^0 ]J \e(fo ) 

Goldman theory [19] results in the following dispersion 
relation 

{e(k, ijj) + xo,o) (^(^ - 2/co, ^ - 2cJo) + X-2-2) = Xo -2^ 

(2) 

where 

^ _ 2^ 4:N'^Jn{zi)Jn{zs) 

""^n, ^^ziZs{{uj-kvo)'-N^ujl)' 

e{k^uj) is the linear plasma permittivity, Jn are the 
Bessel functions of arguments Zj = {k -\- jko)a and ut is 
the space- aver aged density of trapped particles. In order 
to calculate the growth-rate of the mode 4-0 as the func- 
tion of the bounce frequency of BGK-wave, we should 
take into account the contribution of only one resonance 
N = 1. Numerical solution of such a dispersion relation 
is presented in Fig. [4| From this flgure we notice that 
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Figure 4: Dependence of the growth-rate of sideband insta- 
bility for the mode 4-0 on the bounce frequency of trapped 
electrons in a one-dimensional BGK-wave. 

for the observed amplitude of the BGK-wave the best 
agreement between theoretical and simulation results is 
achieved, if the effective density of trapped particles is 
slightly lower than the beam density: Ut = 0.9n5. If for 
the same amplitude of nonlinear wave we take into ac- 
count the contribution of the resonance A/" = 3 in the dis- 
persion equation, we are also able to predict the growth- 
rate of the mode 5-0, which is found to reach the value 
T/ujp 0.015. Besides longitudinal modes 4-0 and 5- 
0, which can be predicted theoretically, we also observe 
instabilities of oblique modes 4-m and 5-m with small 
propagation angles. 

As to the nonlinear saturation of sideband instabilities, 
in the time interval ujpt = 395.6 ^ 659.4 (Fig. [3| this 
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Figure 5: (Color online) Evolution of the beam density map (top row) and beam phase spaces {x,px) (middle row) and {x,py) 
(bottom row) in the case Q = 0.5. Both phase spaces are obtained by averaging on y inside the region restricted by horizontal 
lines in the beam density map. Lines in the phase space {x,px) show px corresponding to phase velocities of different modes: 
3-0 (sohd), 4-0 (dashed) and 5-0 (dash-dotted). 



process is accompanied by transformation of the primary 
BGK-wave with n = 3 to the other nonlinear wave with 
n = 4. From the motion of beam particles in the phase 
space (x^Px) we conclude that the wave energy comes to 
the secondary mode 4-0 not only from the kinetic energy 
of the beam that is decelerated by this slower wave, but 
also from the primary unstable mode 3-0 that loses most 
of its energy during nonlinear wave- wave interactions. In 
the later stage, plasma nonlinearities come into force and 
the BGK-wave with n = 4 is dissipated due to the local 
(in y) buildup of the modulation instability. 



B. Weak magnetic field ft = 0.6 

In the weak magnetic field, due to decreasing in 
growth-rates of oblique instabilities, the main role at the 
nonlinear stage of beam-plasma interaction is played by 
modes 3-0 and 3-1. The exponential growth of these 
modes is saturated by beam trapping. In the chosen local 
region, this nonlinear process can be considered as one- 
dimensional. Indeed, Fig. [5] shows that at the moments 



ujpt = 115.6 and ujpt = 252.5 beam particles mix inside 
the trapping phase-space region of the 3-0 mode without 
a visible increase in their transverse momentum. 

Due to further damping of the primary mode 3-1, we 
observe the formation of almost one-dimensional BGK- 
wave. Since the amplitude of this nonlinear wave is com- 
parable with that in the previous case, secondary insta- 
bilities of almost longitudinal sideband modes with n = 4 
and n = 5 grow with the same growth-rates as in the 
isotropic plasma [Fig. [2|b)]. Nonlinear saturation of 
these modes, however, is not accompanied by the com- 
plete wave energy transfer from primary unstable modes 
to sideband modes. As one can see from Fig. [5] demon- 
strating snapshots of the phase space {x^Px) at the mo- 
ments Upt = 422 and Upt = 688.3, beam particles move 
simultaneously in fields of both primary (n = 3) and sec- 
ondary modes (n = 4). The beam spread in transverse 
momentum, increasing due to excitation of oblique side- 
band modes 4-1 and 4-2, is found to be much smaller 
than in the case of isotropic plasma. 

Thus, in the weakly magnetized plasma the nonlinear 
stage of the beam-plasma instability starts with almost 
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Figure 6: (Color online) Evolution of the Ex map (top row), the beam density map (middle row) and the phase space of the 
beam (x,px) (bottom row) in the case Q = 2. The phase space is obtained by averaging on y inside the region restricted by 
horizontal lines in the beam density map. Lines in the phase space {x,px) show px corresponding to phase velocities of different 
modes: 3-0 (solid), 3-1 (dashed), 3-2 (dash-dotted) and 3-3 (dotted). 



one-dimensional beam trapping. The resulting BGK- 
wave propagating parallel to the magnetic field drives 
sideband instabilities of various upper- hybrid modes. In 
contrast to the previous case, nonlinear interactions be- 
tween primary and secondary modes do not result in the 
substantial energy transfer and dominant modes 3-0 and 
4-0 have comparable amplitudes in the quasi-stationary 
state. 



C. Strong magnetic field ft — 2 

In the strong magnetic field the unstable spectrum nar- 
rows significantly. Beam trapping continues to be the 
main nonlinear process resulting in saturation of domi- 
nant unstable modes 3-0 and 3-1. The contribution of 
oblique mode becomes herewith more pronounced. In 
Fig. [6] demonstrating snapshots of the electric field map 
Ex{x,y), the beam density map ni}{x^y) and the beam 
phase space {x^Px)^ one can see that trapping of beam 
particles in different regions along y starts at different 
moments of time. By the moment ujpt = 187.1, the non- 
linear BGK-wave has not yet been formed, but primar- 



ily stable modes 3-2, 3-3 and 4-0 have already grown 
up to large amplitudes. All of these modes get into the 
Cherenkov resonance with beam particles, but their rapid 
intermittent growth does not allow to clarify unambigu- 
ously what excitation mechanism (Cherenkov or nonlin- 
ear) plays the main role. 

Fig. |2|c) shows that, upon saturation of secondary 
unstable modes, the wave spectrum is dominated by 3- 
m modes with the prevailing role played by the oblique 
mode 3-2. As to determine the branch of plasma oscil- 
lations, to which secondary modes with m > 2 belong, 
we should find real parts of their frequencies. This can 
be done by using the relationship between different field 
components E^^ Ey and B^: 

'^(k) = {k,Ey{k) - kyE^ik)) . 

Substituting simulation results into this formula, we iden- 
tify secondary modes as electrostatic whistler waves. 

Dominance of plasma modes with n = 3 is also visible 
in the phase space of the beam. In Fig. [6] at the mo- 
ment ujpt = 252.5, typical vortex structures arising from 
trapping of beam particles by different 3-m modes are 
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observed clearly. Due to the difference in parallel phase 
velocity, trapping regions of these modes move relative 
to each other, but when bottoms of potential wells have 
the same x-coordinates, electric fields of excited modes 
are summed coherently and form the large potential well 
resulting in phase space mixing of beam particles from 
different trapping regions. From spatial structures of 
the electric field and the beam density at the moments 
oOpt = 252.5 and Upt = 715.9 we see that the spatial re- 
gion of coherent wave- wave interaction is localized in the 
transverse direction and moves along y. Intense phase- 
space mixing of beam particles inside this region makes 
the nonlinear interaction between plasma modes with the 
same parallel wavenumbers k\\ to be more efficient than 
interaction between modes with different n. 

Thus, in strongly magnetized plasmas, nonlinear wave- 
wave interaction via shared trapped particles is marked 
by the prevailing energy transfer from linearly unstable 
modes to initially stable oblique whistlers with the same 
parallel wave number. In spite of excitation of oblique 
plasma waves, the beam spread in transverse momentum 
remains unchanged for the whole period of beam-plasma 
interaction. 



V. DISCUSSION OF SIMULATION RESULTS 

Detailed two-dimensional PIC simulations show that 
collective interaction of a cold low-density electron beam 
with an isotropic plasma evolves according to the follow- 
ing scenario. Prevailing growth of oblique linear instabil- 
ities is followed by trapping of beam particles in oblique 
directions and substantial spreading of the beam in trans- 
verse momentum. The linear growth-rate of the filamen- 
tation instability in our case is much smaller than growth- 
rates of oblique and two-stream instabilities, but from 
the certain time filamentary perturbations demonstrate 
rapid growth due to nonlinear interaction of primary 
unstable modes. Simultaneous buildup of two-stream, 
oblique and filamentation instabilities results in the sit- 
uation when the beam appears to be divided into sepa- 
rate bunches localized spatially in both longitudinal and 
transverse directions. Current layers constructed from 
these bunches merge into one thus forming the nonlinear 
BGK-wave localized in y. This nonlinear equilibrium is 
found to be unstable against oscillations with frequencies, 
which in the reference frame of BGK-wave equal to the 
harmonics of the bounce frequency of trapped electrons. 
Growth-rates of secondary sideband modes propagating 
parallel to the beam direction are well predicted by the 
one-dimensional theory. Saturation of sideband instabil- 
ities is characterized by almost complete energy trans- 
fer from primary unstable modes to almost longitudinal 
Langmuir modes with shorter wavelengths. It results in 
the formation of the BGK-wave with /cy > ujp/vb^ which 
is followed by dissipation of its energy due to the modu- 
lation instability. 

In the weak magnetic field, oblique instabilities cease to 



dominate in the linearly unstable spectrum and the non- 
linear stage of beam-plasma interaction starts with al- 
most one-dimensional beam trapping followed by the for- 
mation of one-dimensional BGK-wave. In contrast to the 
case of isotropic plasma, the beam does not acquire sig- 
nificant spread in transverse momentum at this stage. Al- 
though sideband instabilities of short-wavelength modes 
4-m and 5-m grow with the same rates as in the isotropic 
plasma, saturation levels of these secondary upper-hybrid 
modes appear to be lower. In this case, the prevailing 
energy transfer at the saturation stage remains to be di- 
rected along but nonlinear interaction between pri- 
mary and secondary modes becomes less efficient. 

The strong magnetic field does not affect essentially 
the initial stage of beam trapping by linearly unstable 
modes with small propagation angles, but modifies sub- 
stantially the spectrum of secondary modes arising due 
to nonlinear wave-wave interaction via the beam nonlin- 
earity. In this case, the wave energy of primary unsta- 
ble modes is transferred to oblique modes with the same 
parallel wavenumber. It means that the energy flow is 
directed predominantly along k_\_. 

Q=0.5 




Figure 7: (Color online) Phase velocities of linearly stable 
(blue) and linearly unstable (red) modes in the beam-plasma 
system in terms of momentum of resonant electrons. Bold 
arrows indicate the prevailing direction of energy flow. 

In order to understand why the direction of spectral 
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energy transfer at the nonlinear stage of beam-plasma in- 
teraction changes drastically with the increase of external 
magnetic field, let us consider some dispersion character- 
istics of plasma modes that get into the Cherenkov res- 
onance with the trapped electron beam. Let us express 
phase velocities of these modes in terms of momentum of 
resonant electrons 

rrieC l^V^k/ J 

Fig. [7| presents the dependence of this momentum on the 
transverse wavenumber for different branches of plasma 
waves and different magnetic fields. Squares show posi- 
tions of modes that can be resolved in our periodic system 
and colors indicate amplitudes of these modes at the end 
of simulation run. The main difference between cases 
Q = 0.5 and 1] = 2 is that in the weak magnetic field 
the unstable beam branch couples to the upper-hybrid 
branch whereas in the strong magnetic field it couples 
to the whistler branch. In the former case the phase ve- 
locity of unstable plasma modes with the same parallel 
wavenumber increases with the increase of k_\_, whereas 
in the latter case it decreases even after the transition 



from unstable beam modes to stable whistlers. Thus, the 
observed distribution of wave energy between different 
plasma modes can be explained, if nonlinear wave-wave 
interactions due to the beam nonlinearity are governed by 
the following rules: (i) interaction between modes is more 
efficient, if these modes have the same k\\; and (ii) wave 
energy is transferred to modes with lower phase veloci- 
ties. It explains why, in the weak magnetic field, oblique 
unstable modes 3-m, despite of their large growth-rates, 
are saturated at rather low levels and why the mode 3-1 
loses its energy at the later nonlinear stage. In the ab- 
sence of slower modes with /cy = ojp/vb^ wave energy here 
is transferred to the modes with k\\ > uOp/v}) arising due 
to sideband instabilities. In the strong magnetic field, the 
efficient energy transfer from primary unstable modes 3-0 
and 3-1 to other modes with the same k\\ and lower phase 
velocities becomes possible. That is why the spectrum of 
plasma oscillations in the relaxed beam-plasma system 
in this case is dominated by oblique whistler modes. 

This work is supported by President grant NSh- 
5118.2012.2, grant 11.G34.31.0033 of the Russian Feder- 
ation Government and RFBR grants 11-02-00563-a, 11- 
01-00249-a. 



[1] A. Bret, L.Gremillet, D. Benisti and E. Lefebvre, Phys. 

Rev. Lett. 100, 205008 (2008). 
[2] A. Bret, L. Gremillet and D. Benisti, Phys. Rev. E 81, 

036402 (2010). 

[3] L. Gremillet, D. Benisti, E. Lefebvre, A. Bret, Phys. Plas- 
mas 14, 040704 (2007). 

[4] J. T. Frederiksen and M. E. Dieckmann, Phys. Plasmas 
15, 094503 (2008). 

[5] X. Kong, J. Park, C. Ren, Z.M. Sheng, and J.W. Tonge, 
Phys. Plasmas 16, 032107 (2009). 

[6] A. Karmakar, N. Kumar A. Pukhov, O. Polomarov and 
G. Shvets, Phys. Plasmas 15, 120702 (2008). 

[7] A. Karmakar, N. Kumar A. Pukhov, O. Polomarov and 
G. Shvets, Phys. Rev. E 80, 016401 (2009). 

[8] A. Bret, L.Gremillet and M.E. Dieckmann, Phys. Plas- 
mas 17, 120501 (2010). 

[9] M. Oppenheim, D.L. Newman, and M.V. Goldman, 
Phys. Rev. Lett. 83, 2344 (1999). 
[10] L. Muschietti, I. Roth, C.W. Carlson, and R.E. Ergun, 

Phys. Rev. Lett. 85, 94 (2000). 
[11] D.L. Newman, M.V. Goldman, M. Spector, and F. Perez, 



Phys. Rev. Lett. 86, 1239 (2001). 
[12] N. Singh, S. M. Loo, and B. E. Wells, J. Geophys. Res. 

106, 21183, DOL 10.1029/2001JA900056 (2001). 
[13] T. Umeda, Y. Omura, T. Miyake, H. Matsumoto, and 

M. Ashour-Abdalla, J. Geophys. Res. Ill, A10206, DOL 

10.1029/2006JA011762 (2006). 
[14] T. Umeda, Phys. Plasmas 15, 064502 (2008). 
[15] LV. Timofeev and A.V.Terekhov, Phys. Plasmas 17, 

083111 (2010). 

[16] A.V. Burdakov, A. A. Ivanov, E.P. Kruglyakov, Plasma 

Phys. Control. Fusion 52, 124026 (2010). 
[17] A.V. Burdakov, A.V. Arzhannikov, V.T. Astrelin, 

A.D. Beklemishev, A. A. Ivanov, LA. Kotelnikov, E.P. 

Kruglyakov, S.V. Polosatkin, V.V. Postupaev, S.L. Sinit- 

sky, LV. Timofeev, and V.P.Zhukov, Fusion Sci. Technol. 

59(1T), 9 (2011). 
[18] W.L. Kruer, J.M. Dawson, and R. Sudan, Phys. Rev. 

Lett. 23, 838 (1969). 
[19] M.V. Goldman, Phys. Fluids 13, 1281 (1970). 



